Frustration effects in rapidly rotating square and triangular optical lattices 
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We discuss the ground state of the two-dimensional Bose- Hubbard (BH) Hamiltonian, relevant for 
rotating gaseous Bose-Einstein condensates, by employing U (1) quantum rotor approach and the 
topologically constrained path integral that includes a summation over U (1) topological charge. We 
derive an effective quantum action for the BH model, which enables a non-perturbative treatment of 
the zero-temperature phase transition. We calculate the ground-state phase diagram, analytically 
deriving maximum repulsive energy for several rational values of the frustration rotation parameter 
/ = 0, 1/2, 1/3, 1/4, and 1/6 for the square and triangular lattice, which improves upon previous 
theoretical treatments. The ground state of the rotating Bose-Einstein condensates on a triangular 
lattice appears to be most stable against the effects of rotation. Performed calculations revealed 
strong dependence of the critical ratio of the kinetic energy to the repulsive on-site energy, that 
separates the global coherent from the insulating state, on topology of the lattice. 

PACS numbers: 05.30.Jp, 03.75.Lm, 03.75. Nt 



I. INTRODUCTION 

The merging of atomic and condensed matter physics 
since the experimental realization of Bose-Einstein 
condensation""^ has opened exciting new perspectives for 
the creation of novel quantum states. Especially, systems 
of ultra-cold atoms confined in optical lattices^El facili- 
tate an experimental environment, where a rich variety 
of quantum many-body models can be implemented in a 
wide range of spatial dimensions, geometries, and particle 
interactions. Surprisingly, the quantum phase transitions 
in systems under uniform magnetic field can be also ana- 
lyzed considering rotating Bose-Einstein condensates^"21 
trapped in a two-dimensional (2£') lattice potential. In 
a frame of reference rotating about the z-axis with angu- 
lar velocity Q, the kinetic term in Hamiltonian is equiva- 
lent to that of a particle of charge q experiencing a mag- 
netic field B with qB = 2mU, where m is the mass of 
the particle.EK^ This connection shows that the Coriolis 
force in the rotating frame plays the same role as the 
Lorentz f orce o n a charged particle in an uniform mag- 
netic field .U^JJ^ The presence of angular velocity induces 
vortices in the system described by the rotation frustra- 
tion parameter / (= ma^n/TrH, with a being the lattice 
spacing) . The parameter / can be also expressed in terms 
of the recoil energy En a,s f — TThil/2Efi. Of special in- 
terest are cases when / = p/q, with p and q being the 
rational numbers. Frustration occurs in this system be- 
cause two different area scales are in competition. One 
characteristic area is the unit cell a of considered lattice. 
The other ■nh/mVt is associated with the rotation of the 
lattice. We can use the notion "magnetic field" and "ro- 
tation" interchangeably, assuming that a harmonic con- 
finement potential is applied to cancel the centrifugal ef- 
fects of rotation. Therefore, the nexus of condensed mat- 
ter and optical physics is transparent (since the effects 



of magnetic field/rotation have the same mathematical 
structure) and different systems can mimic each other. 

The progress in setups contrivance used for creation 
of a rotating optical lattice led to systems with different 
geometries like square (D) or triangular (a) that can be 
analyzed in the strongly interacting regime.— Up to now 
two experimental strategies have been developed.!^ The 
first one is based on direct imrarinting a phase shift on 
the macroscopic wave function.ES The second approach 
called "stirring" is an adaptation of the rotating bucket 
experiment to a gas of trapped bosons.C'*^^ The lat- 
ter method fails when angular velocity is comparable to 
trapping frequency D. ~ uj_\_. However, several groups 
have found a way to circumvent the problem of the cen- 
ter of mass expulsion occurring at fi ^ u)_\_ and one can 
achieve i7 — 1.05ajj^.'^ Mott-insulator (MI) - superfiuid 
(SF) transition boundary obtained by using a Gutzwiller- 
type variational wave function revealed the complexity 
of the dependence of phase boundary on the effective 
magnetic field/rotation, refiecting the self-similar prop- 
erties of the single particle energy spectrum. --Mean field 
theory calculations determined that the linear eigenvalue 
equation characterizing the Mott lobe also characterizes 
the Hofstadter butterfiy spectrum. From this authors 
determined an expression for the Mott-lobe boundary.^iSl 
Despite the several theoretical approaches to the problem 
of strongly interacting bosons in rotating lattices many 
questions still remain open and unsolved. 

The aim of this work is to study the superfiuid to Mott- 
insulator zero-temperature phase transition by means of 
the Bose-Hubbard (BH) model in two-dimensional rotat- 
ing optical condensates with different geometries. We ad- 
dress the question of evolution of the ground state phase 
diagram for the system with various angular velocities 
for square and triangular lattices. The inherent difficulty 
of dealing with BH Hamiltonian appropriate for strongly 



correlated bosons originates from the non-perturbative 
nature of the model and the presence of rotation. To elu- 
cidate the quantum phase transition in optical lattices, 
where the kinetic energy scale is less than the dominat- 
ing interaction energy and angular velocity is compara- 
ble to the recoil energy, we have adopted a theoretical 
approach for strongly interacting fermions^"- to the BH 
model in a way to include the effects of particle num- 
ber fluctuations and make the qualitative phase diagrams 
more quantitativepl] To facilitate this task, we employ 
a functional integral formulation of the theory that en- 
ables to perform functional integration over fields defined 
on different topologically equivalent classes of the U (1) 
group, i.e., with different winding numbers. An inclusion 
of the winding numbers is unavoidable in order to obtain 
a proper phase diagram. The quantum rotor representa- 
tion method we use is deeply rooted in the gauge sym- 
metries of the model. We construct an invariant theory 
introducing an appropriate U(l) gauge transformation. 
The outline of the paper is as follows In Sec. II we 
introduce the model Hamiltonian and the effects of ro- 
tation are discussed in Sec. III. Next, we derive an ef- 
fective U (1) action in the quantum rotor representation 
described in Sec. IV- VII. The aim of Sec. VIII is the 
presentation of the resulting phase diagrams for two- 
dimensional square and triangular Bose-Hubbard sys- 
tems in rotating frame and comparison of our results with 
several numerical and analytical calculations. Finally, 
Sec. IX summarizes our results. In the Appendix, we 
give an analytical derivation of density of states (DOS) 
in closed form for several rational values of / = p/q in 
D and A lattice. Moreover, the connection between the 
DOS and Hofstadter butterfly is shown. 



II. MODEL 

In optical lattices the two main energy scales are set 
by the hopping amplitude proportional to t (that sets 
the kinetic energy scale for bosons) due to the parti- 
cles tunneling, and the on-site interaction [/ > 0. For 
t > U the phases of the superfluid order parameter on 
individual lattice sites are well deflned. On the other 
hand, for sufficiently large repulsive energy [/, the quan- 
tum phase fluctuations lead to complete suppression of 
the long-range phase coherence even at zero tempera- 
ture. The competition between the kinetic energy, which 
is gained by delocalizing bosons over lattice sites and the 
repulsive interaction energy, which disfavors having more 
than one particle at any given site, can be modeled by 
the following quantum Bose-Hubbard HamiltoniarP^ 



relations [a^, a , 



alar is the boson num- 



H = — - 2_^nr (n-7- — 1) — 2^ trr'al.ar' — fJ,2_^^r, (1) 



where al. and a^' stand for the bosonic creation and anni- 
hilation operators that obey the canonical commutation 



ber operator on the site r, and the chemical potential 
/i controls the number of bosons. Here, {r,r') identifies 
summation over the nearest-neighbor sites. Furthermore, 
trr' is the hopping matrix element with dispersion t, ' . 



III. EFFECTS OF ROTATION 

In the fast rotation regime the physics of Bose-Einstein 
condensates is very reminiscent of that of charged par- 
ticle in magnetic field. If the centrifugal term can be 
compensated by the trapping frequency in the plane per- 
pendicular to the rotation axis {Q = uj±), so only the 
Coriolis term is left, we have situation which formally is 
equivalent to the Lorentz force exerted by uniform mag- 
netic field on charged particle. Experimentally a region 
of fast rotations up to 17 = 1.05uj±^ that suits our the- 
oretical predictions can be achieved. 

An angular velocity enters the Hamiltonian Eq. (fl]) 
through the Peierls phase factor according to 



t^ 



t^r' CXp 



2m 



A-dl 



(2) 



where A{r) — H x r is the equivalent of a magnetic 
vector appears from the rotation and n = h/m is the 
quantum circulation unit. Thus, the phase shift on each 
site is determined by the vector potential A (r) and in 
typical experimental situations can be entirely ascribed 
to the external magnetic field/angular velocity. We as- 
sume throughout this paper that the model in Eq. (fl]) is 
defined on a lattice with lattice spacing a = 1. From Eq. 
(I2]), it follows that the properties of the system will be 
periodic with a period corresponding to 



/ = 7rhn/2ER 



(3) 



per plaquette. Of special interest are the values of the 
angular momentum which correspond to rational num- 
bers oi f — 1/2,1/3,1/4,... Since all properties of the 
Hamiltonian Eq. (fTj) are invariant under f ^> —f and 
also under / — > / -I- 1, it is sufficient to consider / in the 
range < / < 1/2 that can be reached experimentally. 



IV. DESCRIPTION OF THE METHOD 

We write the partition function of the system switching 
from the particle-number representation to the conjugate 
phase representation of the bosonic degrees of freedom 
using the bosonic path-integral over the complex fields 
Ur (t) depending on the "imaginary time" < r < /3 = 
l/k^T with T being the temperature 



[VaVa] exp 



/9 



drU (r) 



E 



d 



drar (r) — a^ (r) 



(4) 



We decouple the interaction term in Eq. (fTl) by a Gaus- 
sian integration over the auxihary scalar potential fields 



Vr (r) = T4o + K' {r) , 



with static 



Vr0^r^Vr{uJu=Q) 



and periodic part 



+00 



c.c, 



(5) 



(6) 



(7) 



where tUy = 2'kv/I3 {v = 0,±1,±2, ...) is the Bose- 
Matsubara frequency. We observe now that the BH 
Hamiltonian has a local U (1) gauge symmetry, when 
expressed in terms of the underlying boson variables. 
This points out a possibility of an emergent dynami- 
cal U (1) gauge field as a fiuctuating complex field at- 
tached to bosonic variables, which is dynamically gen- 
erated, by interacting bosons. Thus, the periodic part 
V^ (t) = V^ (r -I- /3) couples to the local particle num- 
ber through the Josephson-like relation 0^ {t) = V^ (r) , 
where 



(t>r{T) 






-•t>r.{T) 



ld_ 
i dr 



AAr) 



(8) 



The quantity (p (r) is the U (1) phase field and satis- 
fies the periodicity condition 0^ (/3) = ^^ (0) as a con- 
sequence of the periodic properties of the V^ (t) field in 
Eq. 0. 



V. PARAMETRIZATION OF THE BOSON 
FIELD AND THE ORDER PARAMETER 

We perform the local gauge transformation to the new 
bosonic variables 



Or (t) 




■ gi0,(r) 







\ br (r) 


Or (t) 




t 


;-»0r.(r) 




brir) 



(9) 



that removes the imaginary term —i L drcjir (t) n^ (t) 
from all the Fourier modes except at zero frequency. 
From the above we deduce bosons have a composite 
nature made of bosonic part b^ (t) and attached "fiux" 
exp[i(j)r (t)]. We parametrize the boson fields br (t) — 
bo + b^ (r) and incorporate fully our calculations to the 
phase fiuctuations governed by the gauge group U(l). 
Assuming nonfluctuating amplitude at low temperatures 
br (t) = &01 we drop the corrections, which was proved 
to be j ustified in the large U/t limit we are interested 
jj^| 2i | 23 | rpjjg calculation of 60 is postponed to the next 
section. 

It is very convenient to define the order parameter 




Figure 1: (Color online) The phase variable's part tpB of the 
order parameter \['s Eq. ( 10 1 for rotating triangular lattice 
with the rotation frustration parameter / = 1/4. Below the 
surface the phase coherent state (SF) takes place. Flat region 
means the incompressible Mott-insulator (MI) phase. The up- 
per panel is the density plot of the surface in order to highlight 
the interaction t/U- chemical potential jj,/U dependence. 



which signals the emergence of the superfluid phase and 
vanishes in the Mott-insulator state. The SF state is 
characterized by spontaneously breaking of the U (1) 
symmetry of the Bose-Hubbard Hamiltonian. Note, that 



a nonzero value of the amplitude 5o in Eq. ( 10 ) is not suf- 



^i 



(ur (r)) = (br (r) exp [if/)^ (r)]) = 6o^i 



(10) 



ficient for superfiuidity. To achieve this, also the phase 
variables in Eq. (10 1, must become stiff and coherent, 
which implies ipB 7^ (see Fig. IT] and Fig. l2|. In the 
symmetry breaking state, with a finite expectation value 
of ar (t), different phases (j)i,(j)2, ■■■,(t>n of the condensate 
lead to degenerate ground states (Fig. 2|. If we change 
the phase of the condensate in a large but finite region, 
then, locally, the system is still in one of the degenerate 
ground states. Slow changes of the phase result in the ap- 
pearance of the low energy excitations that correspond 
to fiuctuations among the degenerate states. Based on 
this picture we identify &' as high energy fast fiuctua- 
tions (and drop them in calculations) contrary to the low 




Figure 2: (Color online) The order parameter ^s can achieve 
nonzero value when both amplitude bo and phase (j> are posi- 
tive. The fluctuations of the amplitude b' in the Matsubara 
time in low-temperature limit T — >■ are dropped in our ap- 
proach. 



energy fluctuations described by 



VI. PHASE ONLY ACTION 

By integrating out the auxiliary static field Vro we cal- 
culate the partition function with an effective action ex- 
pressed in the form of the propagator G 



(11) 
where fi/U — fi/U + 1/2 is the shifted reduced chemical 

potential. In the above exp I — TrlnG^^ ) = detG and 

the determinant takes the form 



detG 



/ [Dim] exp J - V / 



dr 



Xbr 



dr 



+ fi]S, 



The inverse of the propagator becomes 

K ^G^Ul- KGo 



G — 1 /^~1 



(12) 
(13) 



The explicit value of the amplitude bo in Eq. ( 10 1 



can be obtained from minimization of the Hamiltonian 
dH {be) /dbf) = 0. Therefore, we write 



Go 



hi 



U 



(14) 



K 



e"l'-^^h>e-"^^'^^\ 



(15) 



Expanding the trace of the logarithm we have 



Tr In G"i = -Tr In Go - Tr 



1, 



(kGo) 



-Tr 



{kg)' 



(16) 



partition function Eq. (11) becomes 



with Go and K given by Eq. ( |14[ ) and ([15 1. Finally the 

(17) 



with an effective action expressed only in the phase fields 
variable 



5phasc [0] = / d,T <y^ 







2^«W + 7l*-M 



J2 e'-^-^^^J^^'e-^-^-'^^) 



> , 



(18) 



where the phase stiffness coefficient is given by J^^' ~ 
b^t^^'. The total time derivative Berry phase imaginary 
term in Eq. (18) is nonzero due to topological phase 
field configurations with (pj. (/3) — (pr (0) = 2nmr {rrir — 
0, ±1,±2...) that result in topological ingredients to the 
correlator we will see below. Therefore, we concentrate 
on closed paths in the imaginary time (0, /3) labeled by 
the integer winding numbers m^. The path-integral 



[I?0] 



E 

{rrir-} 



p^m 



4>r{T)-\-2-7TniT. 



[VCP{T)] 



(19) 

includes a summation over rrir and in each topological 
sector the integration goes over the gauge potentials. 

To proceed, we replace the phase degrees of freedom by 
the unimodular scalar complex field ipj. which satisfies the 
quantum periodic boundary condition -0^ (/3) — -07- (0) . 
This can be conveniently done using the Fadeev-Popov 
method with Dirac delta functional resolution of unitjffHj 
where we take iji^. ss & continuous but constrained (on 
the average) variable to have the unimodular value 



1 = / [DipVi)] 5 



I^.(t)| -iv 



i't'T-ir) 



(^r 



-i4>T-('r) 



(20) 



where N is the number of lattice sites. Introducing the 
Lagrange multiplier A, which adds the quadratic terms 
(in the ipr fields) to the action Eq. ( 18 1, we can solve for 



the constraint. The partition function can be rewritten 
to the form 



Z = 



+200 



2Tri 



-Npj^{\) 



(21) 



with the free energy density T 
T = -\ 



\d.Z/ PN given by: 

/3 



— In / [D%l>Vil:] exp < V / drdr 
^P J 1 , ^0 






(22) 



where I^^,' — 1 if r, r' are the nearest neighbors and 
equals zero otherwise, and 



Irr' (t, t') = (exp {-i [(/)r- (t) - 0^' (t')]}) 



(23) 



is the two-point phase correlator associated with the or- 
der parameter field, where (...) denotes averaging with 



respect to the action in Eq. (18 1. The final form of the 



correlator with topological contribution (summation over 
integer winding numbers) 



(r,r') 



Srr' exp I 

E{m,}exp 



X 51 l^^^P 

{rrir-} 



up 

2 

2 



u) 



rrir H 

"^ U 



[-C/(m. + i^)(r-r')]} (24) 



X exp 
after Fourier transform, can be written as 

UI3 



where 



1 4 
Z~oU 



+ 00 

E 



exp 



771 - 



ElV 
u) 



m— — OQ 



1 - 4 (m 



u 



^ifY 



+ 00 



Zo = ^ exp -up (m + ji/Uf /2 



(25) 



(26) 



m— — oo 



is the partition function for the set of quantum rotors. 
The form of Eq. ( 25 I assures the periodicity in the imag- 



inary time with respect to ^/U = ^/U + 1/2 which em- 
phasizes the special role of its integer values. The action 



Eq. (18 I, with the topological contribution Eq. (24 1 



after Fourier transform, is written as 

5eff [V', ■^l = ]^ E "^ki^r^ ^ {i^u) -IpViv, (27) 

where F^^ [uj^) =A — Jk + 7~^ (wiz) is the inverse of the 
propagator. 



VII. CRITICAL LINES 

Within the phase coherent state the order parameter 
■03 is evaluated in the thermodynamic limit A'^ — )■ oo by 
the saddle point method 5F /5\ = and the unimodular 



condition of the U (1) phase variables translates into the 
equation 



1 






iV/3f-A-,/k + 7-M^.) 



(28) 



The phase boundary is determined by the divergence of 
the order parameter susceptibility Fk=o (w^=o) = 



An 



Jp7q--f7"'(c^.=o)-0 



(29) 



which determines the critical value of the Lagrange pa- 
rameter A = Aq, that stays constant in the whole global 
coherent phase. To proceed, it is desirable to introduce 
the density of states for a 2D lattice in the rotating frame 
in the form 



Pp/q 



M)-^i:^(«-^' 



(30) 



n.A 



with tf, ' being the Fourier transform of the hopping ma- 
trix elements. In this context the quantity J™'*^ in Eq. 



p/q 



( 29 1 represents the maximum of the spectrum described 
by the DOS Eq. (30). The problem of computing of 



P /^ [0 reduces effectively to the solution of the Harper 
equation relevant, e.g., to tight binding electrons on a 
two-dimensional lattice with an uniform magnetic flux 
per unit plaquette. In the Appendix, we give an ana- 
lytical derivation of p ;^ (^) in closed form for several 
rational values of p/q. With the help of the above and 
after summation over Bose-Matsubara frequency uj,^, the 
superfluid state order parameter becomes 



1-V^B = 



1 



+ 00 



pS: (0 d^ 



2i{22 



u 



U 



1) ± 

2) U 



{§) 



(31) 

In Eq. ^ V (fi/U) = frac {^l/U) - 1/2, where frac (x) = 

X — [x] is the fractional part of the number and [a;] is the 

floor function which gives the greatest integer less than 

or equal to x: £ = jmax _ c with J™"^ stands for the 
^ ' ^ p/q ^ p/q 

maximum value of the dispersion spectrum t^ ' and z is 

the lattice coordination number. 



VIII. PHASE DIAGRAMS 

The zero-temperature phase diagram of the homoge- 
neous Bose-Hubbard model Eq. (Ill can be calculated 
from Eq. (31) and is shown schematically in Fig. ^ 



as a function of t/U, with the density controlled by a 
chemical potential n/U. At U/t — 7>0, the kinetic energy 
dominates and the ground state is a delocalized super- 
fluid, described by nonzero value of the superfluid order 
parameter ^b t^ 0. At small values of t/U, interac- 
tions dominate and one obtains a seri es of MI lobes with 
flxed integer flUing hb = 1, 2, ...123221 xhe transition be- 



^ 



QMC 
0.05974(3) 



DPT 

0.05909 



MFT 
0.043 



PA 
0.059 



QRA 
0.06719 



Table II: Comparison ol the maximum of the critical value for 
t/U (as a function of the normalized chemical potential ^/C/) 
at the tip of the first {jib = 1) MI lobe for the square lattice 
with several numerical (QMC - quantum Monte-Carlo,'^^ DPT 
- diagrammatic perturbation theory^ and an aly tical works 
(MFT - mean-field theoryP^PA - Pade analysisPS! QRA - our 
calculations using quantum rotor approach). 



tween the SF and MI phases is associated with the loss 
of long-range order. Let us introduce the notation for 
the niaxiniuni of the critical value for parameter t/U (as 
a function of the normalized chemical potential /x/C/) at 
the tip of the first [ub = 1) MI lobe for different lattices 
and frustration parameters / as follows 



□ ,A 



critJ / 



n,A 



(32) 



In Table [III we com pare values of the Xq resu lting 
from several numericaPSESl and analytical studies .'2323 
We found them in good agreement however, mean-field 
theory calculations of the BH model underestimate x^ 
and in the quantum rotor approach there is a slight up- 
ward trend of the boundary towards higher critical values 
of parameter x^ than obtained from numerical calcula- 
tions. The ground state of the rotating Bose-Einstein 
condensates on a triangular lattice appears to be most 
stable against the effect of rotation (see Fig. jsl. The 
stability comes from the higher values of the repulsive 
energy for the triangular lattice. However, if the rotation 
frustration parameter is equal / = 1/3 and 1/2 the ratio 
Xf /a;? of the energy needed to cause loss of the global 
coherent state changes character and is higher for trian- 
gular lattice unlike the cases with / = 0, 1/6 and 1/4 (see 
Fig. Hand Fig. jsl. In the above we choose /i/C/ = 1/2 
because the transition at integer density belongs to the 
universality class of the 2 + 1 dimensional XY model by 
contrast to transition if one cross SF-MI phase boundary 
by variations in the chemical potential. Behavior of the 
maximum repulsive energy a; + in the rotating system 
with / 7^ taken for special value of the ji/U — 1/2 
is non-monotonical in both square and triangular lattice 
(Fig. 6 and Fig. 7|. While critical values for 
ferent lor varioustopologies of the system, t 
seen in the time-of-flight images occurs rather rapidly 
with increasing lattice depth.EI Because the experimental 
parameter Vq/Eh (Vq is the maximum value of the lat- 
tice depth), depends logarithmically on U/t, the small 
changes of the dimensionless depth of the optical lattice 
can cover a wide range of the phase diagram. In order 
to verify the calculated phase boundaries experimentally 
one shall be able to obtain the higher resolution than 
required to distinguish the ub = 1 from ub = 2 tran- 
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Figure 3; (Color online) Phase diagram for square D and 
triangular A lattice (number of particles per lattice site is 
Ub = 1 inside the first and n^ = 2 inside the second lobe 
respectively) with no rotation / = 0. Within the lobes the 
MI phase takes place with ^b = (see also Fig. [I|. 
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Figure 4: (Color online) Phase diagram for square D and 
triangular A lattice (number of particles per lattice site is 
Ub = 1 inside the first and ub = 2 inside the second lobe 
respectively) with rotation frustration parameter / = 1/2. 
Within the lobes the MI phase takes place with 'I's = (see 
also Fig. [T|. 



sition with Xq ' {ub 



2)/xo [ub 



1) w 0.65. With 
increasing number of particles per lattice site (n^ ^'oo), 
the system possesses an exact particle-hole symmetry 
thus, there is no difference between left Jl/U = n — e 
and right /l/f/ = n + e branch of the n-th lobe, where 
e < 0.5. 



IX. CONCLUSIONS 

The physics of strongly correlated bosonic systems is 
the competition between two tendencies of the bosons to 
spread out as a wave and to localize as a particle com- 
bined with a frustration caused by rotation. We pre- 
sented a field-theoretic study of the ground-phase dia- 
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Figure 5: (Color online) The maximum of the critical value 
for t/U parameter (as a function of the normalized chemical 
potential fi/U) at the tip of the first [ub = 1) MI lobe x,. /xf 
for rotating triangular to square lattice. 
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Figure 6: (Color online) The maximum of the critical value 
for t/U parameter (as a function of the normalized chemical 
potential fJ./U) at the tip of the first (ub = 1) MI lobe Xq /xf 
for rotating square lattice. The vertical dashed line marks the 
ratio of the maximum of the critical value for t/U parameter 
for the second to first lobe a;o {ns ~ 2) /xq {jib = 1). 



gram in quantum two-dimensional gaseous Bose-Einstein 
condensates where mentioned emulation takes place. We 
calculated the phase diagram using the quantum rotor 
approach with exactly evaluated density of states for two- 
dimensional lattices with rational magnetic flux/rotation 
frustration parameter / = p/ q for a number of values 
/ = \/q. In systems that are in the global coherent state 
at / = 0, but with the ratio t/U close to the critical value 



{tiu\ 



a rotation can be used to drive the condensates 



into the MI state (Fig. p|. We compare the maximum 
of the critical value for tfU parameter (as a function of 
the normalized chemical potential ///f/) at the tip of the 
first {riB = 1) MI lobe for square lattice with several nu- 
merical and analytical works and found them in a good 



Figure 7: (Color online) The maximum of the critical value 
for t/U parameter (as a function of the normalized chemical 
potential ^/U) at the tip of the first (ns = 1) MI lobe Xq /x ^ 
for rotating triangular lattice. The vertical dashed line marks 
the ratio of the maximum of the critical value for t/U param- 
eter for the second to first lobe Xq (ub = 2) /xq (jib = 1). 



agreement. Note that the dependence of the x ^ ' /xq ' 
from frustration parameter / is non-monotonical (Fig. ^ 
and Fig. I7|. The critical values of the energy needed 
to drive a rotating condensate out of a global coherent 
state change by varying the frustration parameter and 
strongly depends on topology of the lattice. 

The nice feature of presented approach, described in 
details above, is that all the expressions and handling 
are analytic. It is also worth to notice that we pro- 
vide an exact formulas for density of states that can 
be very useful in various situations whenever the mag- 
netic field/rotation is applied to the physical system. To 
our knowledge the analytical expressions for DOS for tri- 
angular lattice were not known in the literature. No- 
tice, we consider only the limit T — 7> since in two- 
dimensional systems with a continuous symmetry the 
long-range order is destroyed by the quantum fluctua- 
tions at finite temperature .l^S! Moreover, we want to em- 
phasize that our approach cannot be used for analysis of 
the Berezinski-Kosterlitz-Thouless transitions since it is 
appropriate only for physical systems where long-range 
order appears. 
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Appendix A: Density of States 



In this appendix we give the explicit formulas for the 



density of states Eq. ( 30 I for square and triangular lat- 
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Figure 8: The lattice we used in our calculations is topo- 
logically equivalent to triangular structure. It appears from 
square lattice when we add bonds between next-nearest neigh- 
bors (n, n) and (n -\- 1, n -I- 1), thus if </> is the flux per plaque- 
tte then for triangular lattice we have 4'/2. 



tice structures with uniform magnetic field/rotation. The 
provided analytical expression can be advantageous in 



evaluating sums over momenta in Eq. (28). Moreover 



the connection between the DOS and Hofstadter butter- 
fly will be shown. 

We start from the dispersion relevant for a lattice that 
can be viewed as more general than square and triangular 
since it includes both cases (Fig. jsl 



ik 



-2t [cos fci -I- a cos fc2 -I- 2/3 cos (fci -I- fc2)] ■ (Al) 



Regarding the above parameters a — and (3 — lead 
to ID chain; a — 1 and /3 = to the square lattice 



tt 



-2i (cos/ci + cosfc2) 



(A2) 



and finally a = 1 and /3 = 0.5 to the triangular structure 



fj. = —2t [cos fci + cos ^2 + cos (fci + fc2) 



(A3) 



The different sign of the f3 parameter changes the parity 
of the DOS that will be very useful in some cases of tri- 
angular lattice in rotating frame. Performing integration 
over momenta 



we get: 



p(a,/?,0 



p (a, /3, = / 2^ K ~ cos fci — a cos fc2 — 2/3 cos (fci -I- fc2)] 

J -TT (27r) 



1 K 



7rV(i+")'-(?-/5)' 



(l+a)"-(g-/3)" 
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l + a + j3>^>l-a- (3 
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l + a- p>i> -l-a + P 
{1+c.f-ii-pf ) f"^"^ l-a-l3>i>-l + a-l3 



i{a+P(.) 



r 



(A4) 



(A5) 



-1/2 



where K (fc) = JJ" dx (l - fc^ sin^ x) is the complete with rotation frustration parameter / = p/q is given by 



elliptic function of the first kind.l^ 



Appendix B: Square lattice in rotating frame 
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p — ifcl 



'''' Mn J 



-t° - 2 cos (fca + 2Tr fn) . 



(Bl) 



(B2) 



The effects of homogeneous magnetic field/rotation 
on particles have many interesting features. While en- 
ergy levels are quantized into Landau levels in 2D uni- 
form space, a very rich structure appears e.g. Hofs- 
tadter butterfiy '^'^ , when the lattice geometry is taken into 
account.!^ If one uses the Landau gauge A (0, x, 0) then 
the dispersion t^ for a square lattice of spacing a = 1 



Equation (Bll is known as Harper's equation and has 



been studied extensively. If integers p and q are chosen 
to represent the angular velocity (with no common factor 
in p and q), then the dependence on the wave vector 
always appears through the generalized structure factor 
In — cosnfci -I- cosnfc2. The density of states given by 



Eq. ( 30 I can be obtained by computing energy bands t^ 
from the eigenvalue equation (Bl I see Table IV 
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- 12ti + 2Ui - 4 - 276^ = 



1. Square lattice without rotation - f — 

In the case of zero rotation the density of states for the 
square 2D lattice reads simply 



Po° (C)=P (1,0,0 • 



(B3) 



Table IV: Energy dispersion for rotating square lattice for 
values / — p/q used in calculations. 



2. Rotating square lattice with / — 1/2 



The calculation of the exact formulae for DOS is 
straightforward, although for large values of q may only 
be done numerically. However, for a number of q values 
of interest it can be calculated analytically with a closed- 

p/q 



form expression for p , (^ as the end result. Below we 



list these cases. 



/.?,= (« = ^p? 



l(e-i) 



e (8 - ^2) . (B4) 



In the above Q (x) is the unit step function. 



3. Rotating square lattice with / — 1/3 



P?/AO = 4^|(^'-2)v/C^ 
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(B5) 
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4. Rotating square lattice with / — 1/4 
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5. Rotating square lattice with / — 1/6 



P?/6(e) = ^ 
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Figure 9: (Color online) a) The density of states for several values of the frustration parameters f — 0, 1/2, 1/3, 1/4, and 1/6 
in square lattice, b) The projection of the DOS on ^ — / surface results in Hofstadter butterfly. Red lines correspond to our 
analytical results. 



v2 sec ( V' ~ TT ) 
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The projection of the analytically calculated density of 
states for several values of rotation on the ^ — / surface 
results in the Hofstadter butterfly (see Fig. l9|. The cal- 
culations of the Hofstadter spectrum feasible for general 
purposes are impractical in applications since the phase 



boundary calculated from Eq. (31 1 strongly depends on 
the structure of DOS. 



Appendix C: Triangular lattice in rotating frame 

For the lattice topologically equivalent to triangular 
(Fig. Is]) the Harper's equation takes form: 
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where n = q/2 for even and n ^ q for odd q. Moreover 



^A 



Mn = -tfc - 2 COS (fc2 + 471 fn) , 



Nn = -e 



— iki 



g27r/(2n+l)gj(fei+fc2) 



(C2) 

(C3) 



and Nn is a complex conjugation. Now the generalized 
structure factor can be written in form 

7_i_,j = — cosnki ~~ cosnk2 =F cos {nki + 71^2) . (C4) 

From the above we can derive the equations for the en- 
ergy dispersion (see Table VI I for several values of the 
frustration rotation parameter. Analytical results some- 
times come at a price of the complexity of solutions and 
that is the case here. Therefore we omit exact results for 
the dispersions t^ and present them only in the simple 
case oi f — 1/4. 



1. Triangular lattice ^vithout rotation - / = 

Cases with / = (1/3) and / = 1/2 (1/6) are relatively 
simple since the density of states can be easily obtained 



11 



V 


Q 


t^^t- 




1 


ifc + 27f = 




2 


ifc - 2j% = 




3 


il - 9tfc - 6 + 27f = 




4 


i^fe - 6 - 27f = 




6 


i| - 9f fc + 6 - 27f^6 = 



by changing /3 — ?► — /3 or equivalently ^ ^ — ^. 



A 



p^(e)^p(i,o.5,o. 



(C5) 



Table VI: Energy dispersion for rotating triangular lattice for 
values / — p/q used in calculations. 
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Figure 10: (Color online) DOS for triangular lattice without 
rotation f — 0. 



2. Rotating triangular lattice with / = 1/2 
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3. Rotating triangular lattice with / = 1/3 
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where 



4. Rotating triangular lattice with / = 1/4 



2/3 



$ ff ) = 2^/'^ £ f£^ — 9) + ff^ — 3) a/£2 _ 12 (C8) "^^^^ energy dispersion calculated from Harper's equa- 



tion (Bl I can be written in the form: 
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Integrating over the wave vectors belonging to the Bril- 
louin zone, we obtain 
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Figure 11: (Color online) DOS for rotating triangular lattice 
/ = l/2. 
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5. Rotating triangular lattice with / = 1/6 
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Figure 12: (Color online) DOS for rotating triangular lattice 
/ = l/3. 
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Figure 13: (Color online) DOS for rotating triangular lattice 
/ = 1/4. 
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Figure 14: (Color online) DOS for rotating triangular lattice 
/ = l/6. 



